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Abstract 


It is shown that detailed and accurate information about the mass spec¬ 
trum of the massive Schwinger model can be obtained using the technique 
of strong-coupling series expansions. Extended strong-coupling series for the 
energy eigenvalues are calculated, and extrapolated to the continuum limit 
by means of integrated differential approximants, which are matched onto a 
weak-coupling expansion. The numerical estimates are compared with exact 
results, and with finite-lattice results calculated for an equivalent lattice spin 
model with long-range interactions. Both the heavy fermion and the light 
fermion limits of the model are explored in some detail. 
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I. INTRODUCTION 


The Schwinger model [|^^, or quantum electrodynamics in two space-time dimensions, 
has been a source of continuing interest over the years. It is a fascinating model in its own 
right, and also exhibits many of the same phenomena as QCD, such as conhnement, chiral 
symmetry breaking with a U(l) axial anomaly, and a topological 0-vacuum H- It is also 
perhaps the simplest non-trivial gauge theory, and this makes it a standard test-bed for the 
trial of new techniques for the study of QCD. 

Our main purpose in this paper is to explore the usefulness of the strong-coupling series 
approach to this model. It is well known that Euclidean Monte Carlo techniques have proved 
difficult and expensive to apply to models with dynamical fermions, because of the infamous 
“minus sign” problem; and thus, it seems worthwhile to ask whether other techniques such as 
strong-coupling expansions can give useful information in such cases. Modern linked-cluster 
techniques allow one to carry these expansions to very high order. We are particularly 
interested to see if the series approach can describe the non-relativistic or heavy fermion 
limit of these models. 

As a first test, we apply the approach to the Schwinger model, using a Hamiltonian 
lattice framework [Q. We concentrate on a calculation of the spectrum of bound states, 
more specifically the lowest two bound states, for the case of the massive Schwinger model. 
It will be shown that the series approach can give quite detailed and accurate information 
on the spectrum. It will be interesting to see if a similar approach is useful for models in 
higher dimensions. 

Hamiltonian strong-coupling series were hrst calculated for the Schwinger model long 
ago by Banks, Kogut and Susskind [0, and were extended by Carroll et al. 0. Since then, 
however, the method has fallen into abeyance, being eclipsed by the power and accuracy of 
the Monte Carlo method. 

Many other approaches have been made to the bound-state spectrum. In the massless 
case, of course, the model is exactly solvable, as shown by Schwinger |^,|^. It is equiva¬ 
lent to a theory of free, massive bosons. For small fermion mass, one may perturb about 
the zero-mass limit, and obtain a low-mass expansion for the spectrum |^. This expan¬ 
sion has recently been carried to second order by Vary et al. 0 and Adam |T^. In the 
large mass or non-relativistic limit the “positronium” bound states can be solved in terms 
of a Schrodinger equation with a linear Coulomb potential |^. For hxed, finite fermion 
mass, one has to resort to numerical techniques. A quite accurate variational calculation in 
the infinite momentum frame was performed by Bergknoff Finite-lattice Hamiltonian 
calculation were performed by Crewther and Hamer |T^ and Irving and Thomas |T^: we 
include some further hnite-lattice calculations in this paper, mainly as a check on the series 
results. Later on Eller, Pauli and Brodsky applied a “discrete light-cone quantization” 
(DLCQ) approach to the problem, and showed that it gave quite good results, not only for 
the lowest state, but for a whole range of higher excited states. Mo and Perry II have used 
a different light-front field theory approach, together with a Tamm-Dancoff approximation, 
which appears to give excellent results. Tomachi and Fujita Q have used a “Bogoliubov 
transformation method”, which works quite well for small fermion masses. 

In Section 2 of the paper the Hamiltonian lattice formulation of the model is reviewed, 
and the known analytic results on the bound-state spectrum are recalled. It is shown that 
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for free boundary conditions the gauge degrees of freedom can be eliminated entirely from 
the lattice model, leading to an equivalent spin lattice model with long-range interactions. 
This is the formulation which we use to carry out the hnite-lattice calculations. 

Section 3 presents the numerical results. Their extrapolation to the continuum limit 
is discussed, and expansions near both the low-mass limit and the non-relativistic limit 
are analyzed. The numerical results are compared with previous approaches. Finally our 
conclusions are summarized in Section 4. 


II. FORMALISM 
A. Continuum Formulation 

The continuum Lagrangian density takes the standard form 

C =+ ^ - g (1) 


where 


F^y = - duAf, ( 2 ) 

and the Lorentz indices p, z/ = 0 or 1. The coupling g in (1+1) dimensions has the dimensions 
of mass. Choosing the timelike axial gauge 

klo = 0 (3) 

the Hamiltonian is found to be 

H = J + igAi)'il: + m'lp'ip + (4) 

where the electric held E has only one component in one spatial dimension: 

E = F^^ = (5) 

The remaining gauge component is not an independent degree of freedom, but can be 
eliminated if desired, using the constraint provided by the equation of motion (Gauss’ law): 

diE = —diA^ = ( 6 ) 

In the massless case, the theory has been solved by Schwinger [m[, and becomes equiv¬ 
alent to a theory of free, massive bosons, with mass 

— = ^ ~ 0.564 (7) 

9 Vtt 

For small electron mass m/g, one can obtain analytic estimates by perturbing about 
the massless limit. Caroll, Kogut, Sinclair and Susskind P| found that the lowest-mass 
(“vector”) state has mass 
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( 8 ) 


Ml 

9 


Tfl TYl 

+ e^(-) + ■ ■ ■ ~ 0.564 + 1.78(-) + 
9 9 


while the ratio of the next-lowest (“scalar”) mass to the vector mass was 


Ml 


2 - + 
9 


9 


( 9 ) 


where 7 ~ 0.5772... is Euler’s constant. 

These results have been extended to second order by Vary, Fields and Pirner 
Adam 


and by 
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= 0.5642 + 1.781(—) + 0.1907(—)2 + ... 
9 9 9 


( 10 ) 


Adam also found 

Ml 


= 2 


TT^e^'’' ,m 


,m. 


(-)^ + ...^ 2 - 1.5368(-)" + 
<oA g g 


( 11 ) 


differing from the result of Carroll et al. by a factor of 2 ^ = 128. 

In the large-mass or non-relativistic limit, the “positronium” bound states are described 
by a Schrodinger equation with linear potential [[TlJl , 


— -I- ) 4/(a;) = E^>ix) 

m 


( 12 ) 


where the non-relativistic energy is 


E = M — 2m 


(13) 


Eq. (0 may be solved in terms of Airy functions, to give the energies of the lowest vector 
and scalar states as 
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0.642(—)^/3 as 
m 

1.473(—)^/3 as 
m 


m/g —>• 00 


m/g —>• 00 


(14) 


B. Lattice Formulation 

The model can now be formulated on a “staggered” spatial lattice |]^. Let the lattice 
spacing be a, and label the sites with an integer n. Dehne a single-component fermion held 
0 (n) at each site n, obeying anti-commutation relations: 

{(t)\n), 0(m)} = 5mn, {0(n), 0(m)} = 0 (15) 
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The gauge field is defined on the links (n, n + 1) connecting each pair of sites, by 

f/(n, n + 1) = (16) 

Then the lattice Hamiltonian equivalent to Eq.(|) is 

N N 2 N 

H = -^^2 + 1) - h-c.] + m^{-l)^(p\n)(p{n) + ^ W (1'^) 

n=l n=l n=l 

where the number of lattice sites N is even, and the correspondence between lattice and 
continuum fields is 


A( M r / V'upper(a:), n even 
0(n)/Va -^ { I ( ( j j 

IV'iowerla:), n odd 


(18) 


■‘upper” and “lower” being the two components of the continuum spinor), and 

^-9{n) —> —A^{x) 


ag 

gL{n) E{x) 

The gamma matrices are represented by 


(19) 


7 = 


0 


0 


7 = 


0 1 
-1 0 


( 20 ) 


As usual, we have chosen a “compact” formulation where the gauge field becomes an angular 
variable on the lattice, and L{n) is the conjugate ‘spin’ variable 


[0{n),L{m)\ = i5nm 


( 21 ) 


so that L{n) has integer eigenvalues L{n) = 0,±1,±2 ,---. As noted by Banks et al. [0, 
this quantization of electric field (or flux) in one dimension also occurs in the continuum 
Schwinger model, due to Gauss’ law. If one takes the naive continuum limit a —0, the 
lattice Hamiltonian (^) reduces to the continuum Hamiltonian (P, as it should. 

Now define the dimensionless operator 

W = -2-H = Wo + xV (22) 

ag 

where 


n n 

V = (j){n + 1) — h.c.] 

n 

2m 1 


(23) 

(24) 

(25) 
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For X -C 1 one can employ strong-coupling perturbation theory on this model, treating 
Wq as the unperturbed Hamiltonian and V as the perturbation, as discussed by Banks et 
al. 0]. In the strong-coupling limit, the unperturbed ground state |0) is the eigenstate with 

L(n) = 0, (j)\n)(j){n) = ^[1 - (-1)^], alln (26) 

whose energy will be normalized to zero, corresponding to a “hlled Dirac sea”. Banks 
et al. 10] have discussed how to use Rayleigh-Schrodinger perturbation theory to generate 
perturbation series in x for the ground state and excited state eigenvalues of this system, 
and the discussion will not be repeated here. We have used more sophisticated linked-cluster 
techniques [|] to generate high-order perturbation series for these eigenvalues, as presented 
below. 

The lattice version of Gauss’ law is then taken as 

L(n) - L{n - 1) = (P\n)(P{n) - ^[1 - (-1)"] (27) 

which means excitations on odd and even sites create =f1 units of flux, corresponding to 
‘electron’ and ‘positron’ excitations, respectively. 


C. Equivalent spin formulation 


The one-component fermion operators can be replaced by Pauli spin operators at each 
site if we employ a Jordan-Wigner transformation 0|, 


giving 


<Pin) = Y[[i(T3il)](^ (n) 


l<n 


= Y[[-i(T3{l)](T^{n) 

l<n 


Wo = Y.L\n) + ^ + N^I2 

n ^ n 

V — + 1) + li.c.] 


(28) 

(29) 

(30) 

(31) 


The strong-coupling ground-state then corresponds to 

L(n) — 0, cr 3 (n) — —( — I)”, all n 
Next, the gauge field ean be eliminated using Gauss’ law: 

L(n) - L(n - 1) = ]^[(T 3 {n) + (-1)”] 
and a residual gauge transformation: 


a (n) 

l<n 




}a (n) 


(32) 


(33) 


(34) 
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provided that we assume free boundaries 

L(0) = L{N) = 0, where N = lattice sites (35) 

[If periodic boundary conditions are assumed, then there is one extra independent gauge de¬ 
gree of freedom left over, corresponding to the “background” electric held [Q]. The resulting 
Hamiltonian is then 


W = Wo+ xV (36) 

where 

u ^ Nu 

ivo = +^+E 

^ n=l ^ n=l 

N-l 

1/ = ^ [cr’''(?7,)cr“(?7, -I- 1) -I- h.c.] (38) 

n=l 


- Y,{(To{m) + (- 1 )' 

/ m=l 


(37) 


All trace of the gauge held has now disappeared, but instead there is a non-local, long-range 
interaction between the spins in the last term of equation (|37|) , which of course corresponds to 
the long-range Coulomb interaction between charges in the original theory. In the continuum 
limit a ^ 0, a; —*• cxo, the interaction V dominates the Hamiltonian so that in leading order 
the system becomes equivalent to a simple XY model with ground-state energy per site 

— —>-as X —>■ oc (39) 

N TT 

The energy gap to the lowest excited state in the sector of “vector” states will correspond 
to the vector mass, so we expect 


2 2 ^ 

cui — cuo = —^Mi = -Miy/x 

ag 9 


(40) 


and similarly in the scalar or ground-state sector, the minimum energy gap 


2 2 

1X2 — UJq = - = —M 2 \/x 

ag 9 


(41) 


Our aim in this paper is to hnd estimates of these masses Mi and M 2 . 

The equivalent spin model has a total of only 2^ possible conhgurations, and lends itself 
readily to hnite-lattice techniques of analysis. As a check on the series results, we have used 
the Lanczos algorithm to obtain exact results for the low-lying eigenvalues of this system 
at hnite x on lattices of up to Y = 22 sites. The ground-state energy uq is easily obtained 
as the lowest eigenvalue in the sector containing the unperturbed ground-state |0); and the 
hrst excited state energy uji is likewise the lowest eigenvalue in the “vector” state sector, 
corresponding in the strong-coupling limit to the state []^ 


| 1 ) 


1 N-l 

-= Y, [cT+(n)cT-(n-l-l) 

VA „_i 


h.c. 


| 0 > 


(42) 


The second excited-state energy 1 x 2 is the lowest of a “band” of excited states in the vacuum 
sector, corresponding in the strong-coupling limit to the state 


| 2 ) 


N-l 




Y, \(T^{n)a (n-I-1)-I-h.c. |0) 


(43) 
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III. RESULTS AND ANALYSIS 


A. Finite-lattice results 

Exact eigenvalues have been calculated for the equivalent spin Hamiltonian (^) using 
the Lanczos technique for various values of m/g and coupling x, on even lattices from 
iV = 4 up to A = 22 sites. No symmetrization of states was employed: since free boundary 
conditions were chosen, the system does not exhibit translational invariance in any case. 
The calculations are carried out in the sector = 0 which has iV!/[(A/2)!]^ states, 

and the ground-state energy ujq and the “vector” excited-state energy oji are the lowest and 
the second lowest eigenvalues in this sector, respectively. The “scalar” excited-state energy 
002 is slightly more tricky to obtain because there are several other states, corresponding to 
the momentum excitations of the “vector” excited state, in this sector giving lower energy 
than the “scalar” excited state. We have got over this problem by moving in stages from the 
strong-coupling limit x = 0 to the desired coupling value, and “tagging” the desired state as 
that which has maximum overlap with its predecessor at each stage, starting from the state 
|2) in (H). 

To make comparison with series data, it is hrst necessary to extrapolate the hnite-lattice 
data to the bulk limit N ^ oo. The convergence of the hnite-lattice sequences was slow. 
Tests showed that the convergence was polynomial in 1/iV: thus for the ground-state energy 

oJo{N) ^ Neo{(X))+ bo +^ + 0{N~^) as N ^ oo (44) 

where eo(cxo) is the bulk energy per site and Bq is the surface energy term; while for the 
energy gaps 


oo,{N)-ooo{N) 


b', + ^ + 0(N 


-3 


) as iV —>• cx) 


(45) 


This is opposed to the exponential convergence ~ exp(—ciV) one normally expects for peri¬ 
odic boundary conditions. There is a simple explanation for this behaviour^]: on the lattice 
with free boundaries, the excitations have hnite momentum 0{7i/N), and thus their ener¬ 
gies include a kinetic energy correction term 0{{7r/N)‘^). As one approaches x —>■ cx), which 
is a critical point of the lattice model, the hnite-lattice corrections become relatively much 
larger, and extrapolation to the bulk limit becomes progressively more difficult. An example 
is shown in Figure 1 for x = 4, m/g = 0. 

Various sequence extrapolation algorithms HISH^ were tried for estimating the bulk 
limit, such as the alternating VBS algorithm, and the Lubkin and Bulirsch-Stoer algorithms, 
but the most reliable and accurate method in this case seemed to be a simple least squares 
polynomial ht in 1/V, with hctitious errors assigned to the data points. Consistency between 
the estimates at diherent orders allows a crude estimate of the likely error in the hnal result. 
This was the technique used to obtain the estimates used in the rest of this paper. 


^We are indebted to Dr. O. Sushkov for this remark. 







B. Series expansions 


Strong-coupling perturbation series have been calculated for the ground-state energy cuq, 
and the energy gaps {ui — cuq) and (a ;2 — ojq), as functions of the coupling x and mass 
parameter jj,. The full series up to order are presented in the Appendix. 

These series can be analyzed in three different regimes, as follows: 


1. Massless limit, mig ^ 0 

For small m/g, the eigenvalues can be expanded as series in the mass parameter /r, e.g. 
for the vector gap 

(cui - cuo) [x, /i] - 2/i = /o (x) -F /i/i (x) + /UV 2 (a:) H- (46) 

The series have been calculated for /j(x) {i = 0,1, •••,6 ) up to order x^® for the vector 
excited state and order x^® for the scalar excited state. Coefficients for the series /j(x) 
(i = 0,1, 2, 3) are listed in Table I. The expected behaviour in the continuum limit for these 
series is 


fi{ 


X 


aix 


( 1 - 0/2 


as X 


CX) 


(47) 


which would lead to a continuum energy gap 


M — 2m 
9 


= ao/2 -0 ai{m/g) + 2 a 2 {m/gY ^ - 


(48) 


These series must now be extrapolated from the strong-coupling limit x = 0 to the con¬ 
tinuum limit X —»■ 00 . We have employed the standard techniques of integrated differential 
approximants and naive Fade approximants foi’ this purpose, combined with a “match¬ 
ing” technique. Some examples are shown in Figs. 2, 3 and 4. 

Figure 2 shows the ground-state energy per site for m/g = 0 as a function oi y = 1/\/x 
(which is the natural variable to use at weak couplings, as shown by equation (|^) and by 
previous weak-coupling analyses p] , pT| ). It can be seen that the approximants converge 
down to about y ~ 0.5, and are in excellent agreement with the hnite-lattice estimates. 
These can easily be extrapolated to the continuum limit x = 0, where 


1^0 _ 

Nx TT 


(49) 


according to equation (|3^). 

A more interesting example is seen in Figure 3, which shows the vector energy gap as a 
function of y. There it can be seen that the series approximants converge down to about 
y ~ 0.7, and then begin to spray outwards; however, the established behaviour is quite 
smooth as a function of y, and a polynomial £t in powers of y to the series data over the 
range 0.7-1.1 gives a “matching” weak-coupling series £t which extrapolates to a value 


Ml — 2m 


9 


0.56(2) 


as 


y 


0 


(50) 
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which compares well with the expected exact value of 0.5642, from equations (H) and (^Ol). 
It can also be seen that the hnite-lattice estimates converge well down to ?/ ~ 0.2, conhrming 
that the extrapolation to weak coupling is valid. A £t to the hnite-lattice data gives an even 
more accurate estimate of the continuum limit 


Ml — 2m 
9 


0.57(1) 


as 


y 


0 


(51) 


Figure 4 is a similar plot for the energy of the scalar excited state. Here it can be 
seen that the series approximants converge down to ?/ ~ 1 , where the function develops 
a pronounced peak or bump. It is not really possible to tell from the diverging series 
approximants whether the function increases, decreases or remains hat below that point, 
but the hnite-lattice data show that it decreases gently towards the exact continuum value 
1.128. A ht to the hnite-lattice data over 0.2-0 .8 in y gives 


M 2 — 2m 


9 


1.14(3) 


as 


y 


0 


(52) 


in quite reasonable agreement with the exact value; but the best one could do with the series 
data would be to estimate a qualitative value. 


of“'“ ~ 1.25(15). 


(53) 


In a similar fashion, estimates have been obtained for the coefficients Oq, Oi, 02 , Os, in 
equation (|48| ) for the vector and scalar masses, extracted from both the series data and the 
hnite-lattice data. These are summarized in Table II, along with the exact results (where 
known). It can be seen that the numerical estimates are in good agreement with the exact 
results for the vector state, and so the estimate for 03 should also be fairly reliable. For 
the scalar state, there is more structure at small y, as seen in Figure 4, and although the 
estimates for oq and Oi agree quite well with the exact results, those for 02 and 03 are 
virtually worthless. 


2. Finite mjg 


At hxed, hnite m/g, the series have been calculated up to order for the ground-state 
energy cuq, and order for the energy gaps {oji — cuq) and (a ;2 — cuo) - these series are 
available on request. The analysis follows very similar lines to those described above. The 
major diherence is that for m/g 7 ^ 0 , the strong-coupling series expansions are in powers 
of rather than x^. The series are therefore longer, but the convergence is not very 
diherent to that at m/g = 0 . 

At large m/g the series coefficients begin to grow very rapidly, like {m/g)'^ at large 
orders n, and it becomes more difficult to maintain good accuracy down to small values of 
y. Figures 5 and 6 illustrate the behaviour of the energy gap at large m/g. It can be seen 
that any structure has moved to smaller values of y, and is somewhat less pronounced than 
at m/g = 0 . 

Our estimates of the energies Ei/g are shown in Table III and Figure 7, along with earlier 


hnite-lattice estimates of Crewther and Hamer [13], and the light-cone estimates of Eller et 
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al. [^] and Mo and Perry |^. It can be seen that all the estimates agree with each other, 


within errors, except that the results of Eller et al. |]^] run a little too high at small m/g. 
For the vector state, the data match on beautifully to the asymptotic expansions at both 
ends. For the scalar state, the data indicate a peak in the energy at about m/g ~ 0.5. 

Comparing the different estimates, it is noticeable that our current hnite-lattice estimates 
are in fact less accurate than the old estimates of Crewther and Hamer or the even more 
accurate results of Irving and Thomas |^, although there is good consistency between 
them. This is because of the large hnite-size corrections associated with the free boundary 
conditions. More importantly, the series estimates are seen to be in good agreement with the 
current hnite-lattice estimates, and for the vector state they are generally almost as accurate 
(approximately 10%), although for the scalar state they are worse (approximately 15-20%). 
The results of Mo and Perry seem very reliable and accurate, but unfortunately there 
is no assessment of the likely errors in their results. 


3. Non-relativistic limit, m/gf —> oo 


At large m/g, a natural re-arrangement of the series for the energy gaps was suggested by 
Hamer |]T^. Instead of variables fi and x, one can re-arrange the series in terms of variables 
l//i and u, where 


X 


g, 2mg‘^a^ 2m/g 


X 


3/2 


(54) 


Then one hnds the series can easily be expanded in powers of the inverse mass parameter 
l//i, e.g. for the vector gap 


{uji — uJo)[u, p] — 2/i — fo{u) H— fi{u) H—^/ 2 (w) -!-■■■ 

g g^ 


(55) 


The continuum limit can then be approached by hrst letting m/g (or p) —> cx) in such a way 
that mg"^ (or u) remains hnite, and then letting a —0, or m —> cx). If we assume that each 


separate term on the right-hand side of (|55D gives a hnite contribution to the continuum 
energy gap in this limit, then we require an asymptotic behaviour for the fi{u) given by 


fiiu) 


CiU 


(Z+l)/3 


as u 


oo 


which would lead to a continuum energy gap 


F _ “ Q /^\(4m)/3 
g ^2 V2m/ 


(56) 


(57) 


1=0 


The leading term 0{{g/mY^^) has the correct behaviour as predicted by (|T^. 

The series have been calculated for fi (Z = 0,1, • ■ •, 6) up to order for the vector 
excited state and order for the scalar excited state. Coefficients for the series fi{u) 
{I = 0,1,2,3) are listed in Table IV for the two energy gaps. The hrst hve coefficients of 
/o(m) can be obtained from the results of Carroll et al. and were listed previously by 
Hamer |T^. The series for fo{u) is identical with that obtained from a lattice version of 


11 
































the non-relativistic Schrodinger equation (|^, as shown by Kenway and Hamer [^, which 


conhrms that we should obtain at least the leading term in the non-relativistic limit correctly 
by this procedure. 

The series for fi{u) must now be extrapolated from the strong-coupling limit m = 0 to 
the continuum limit u ^ oo. The same techniques were employed for this purpose as for 
the previous series in x. The only question is what to take as a weak-coupling variable. No 
detailed weak-coupling analysis has been done in the non-relativistic limit; but empirical 
tests on fo{u), making comparison with the exact results, appear to show that the best 
extrapolations are obtained using ^ = 1/m as weak-coupling variable. The results are shown 
in Figures 8 and 9. 

For the vector state, it can be seen that the integrated differential approximants converge 
down to z ~ 0.2, and the function behaves very smoothly in z, so that a very accurate 
extrapolation to the continuum limit is possible, giving 


^vector _ l_gl8(2) (58) 

which compares well with the expected exact value of 1.617(= 2^/^ x 0.642). Note that the 
series results allow a more accurate estimate for this quantity than the hnite-lattice data do. 

The results for the scalar state are shown in Figure 9, where it can be seen that the series 
approximants converge only down to z ~ 0.7, but the function is again quite smooth in z, 
so that a fairly reliable extrapolation may be made to the continuum limit, as conhrmed by 
the finite lattice data. The resulting estimate is 

^scalar ^ 3_73('3) ('gg) 

to be compared with the exact value of 3.71(= 2^/^ x 1.473). 

As soon as one tries to go beyond the leading order, however, some major problems arise. 
An example is shown in Figure 10, which shows fi{u)/v?/^ for the vector state. It can be 
seen that the series approximants and the finite-lattice estimates agree very well down to 
very small values of z^ and indicate a divergent behaviour for this quantity. A Dlog Fade 
analysis of the series indicates that f\{u) behaves asymptotically more like u than 
Similar results are obtained for all the {//(m), / > 0}, for both the vector and scalar states. 

It therefore appears that the assumption (|56D is incorrect, and that the terms in equation 
(1^) cannot be analyzed separately (beyond the leading order, at least). It is very likely that 
the limiting behaviour of the energy gaps is non-uniform, and that by letting hrst mjg —> cxo, 
and then a —> 0, we have taken the limits in the wrong order. This gives an incorrect result 
for the ground-state energy, for instance. It follows that the expansion (|57| ) is also incorrect; 
and in fact one would probably expect higher-order corrections to involve integer powers of 
(g/m), rather than (g/mY^^. It is an interesting puzzle how to obtain useful estimates of 
the higher-order corrections from the series data in these circumstances. At this stage, we 
have no answer to the puzzle. 


IV. CONCLUSIONS 

It has been shown that strong-coupling series expansions can deliver quite detailed and 
accurate information about the mass spectrum of the Schwinger model. Using integrated 
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differential approximants or other methods, the series can be continued or extrapolated 
well into the weak-coupling regime, and then matched onto a weak-coupling form, which 
allows reasonably accurate estimates of the continuum limit. This has been conhrmed by 
comparison with both hnite-lattice data and exact results. At hxed, hnite fermion mass the 
estimates for the lowest, vector excited state energy were accurate to about 10 percent. For 
the next lowest, scalar excited state, the eigenvalue shows more structure at weak coupling, 
particularly for small m/g, and the continuum estimates are more qualitative, at about the 
15-20 percent level. 

Estimates can also be obtained for the expansion coefficients of the continuum energy 
eigenvalues in powers of {m/g) about the zero-mass limit, mig —> 0. These numerical 
estimates agree with the exact results of Carroll et al. P], Vary et al. and Adam [|I^, and 
extend them by one order. The series results for these coefficients have an accuracy equal 
to or better than that of the hnite-lattice results. 

A complementary expansion was attempted about the non-relativistic limit, m/g —> oo. 
Series estimates for the leading-order term of the non-relativistic energy were extremely 
accurate in this limit, with an accuracy of 0.2% for the vector state, but problems arose 
for the higher-order corrections. There appears to be a non-uniform limiting behaviour 
in this case, and the natural structure of the strong-coupling series does not predict the 
correct limiting behaviour. We have not resolved the puzzle of how to analyze the series in 
this situation. It would be very interesting to compare the series data with more detailed 
analytic calculations in the non-relativistic limit - we hope to address this problem in the 
future. 

Finite-lattice techniques can also give accurate information about this model, as pre¬ 
viously demonstrated by Crewther and Hamer 


131 and Irving and Thomas 


Exact 

hnite-lattice eigenvalues have been calculated for the equivalent spin Hamiltonian on lat¬ 
tices of up to V = 22 sites, and then extrapolated to the bulk limit by polynomial hts in 
1/N. Fitting these estimates to a weak-coupling form, estimates of the continuum limit can 
be obtained which are accurate to about 5-10 percent for both the vector and scalar states. 
This is actually worse than the accuracy of the old hnite-lattice data of Crewther and Hamer 

m- 

The factor which hindered the hnite-lattice calculations from giving even more accurate 
estimates was the large size of the hnite-lattice corrections, 0(1/V^), which were associated 
with the free boundary conditions. In retrospect, it might be better to choose periodic 
boundary conditions. This carries the penalty that the gauge held cannot be entirely elim¬ 
inated, and one is left with one extra gauge degree of freedom corresponding to Coleman’s 
background electric held or hux, which would have to be truncated in some fashion. The 
advantage would be that the hnite-lattice corrections should be much smaller, and the hnite- 
lattice sequence should converge more rapidly. With today’s computers, one could probably 
expect to obtain virtually exact eigenvalues on lattices up to some 20 sites, and thus obtain 
a substantial improvement on previous hnite-lattice calculations 


13,14 


In higher dimensions, the hnite-lattice approach is hardly feasible, because of the huge 
proliferation of basis states with lattice size. Monte Carlo techniques are the only other 
option in this direction, and their power is still rather limited for models with dynamical 
fermions. It remains to be seen whether modern series expansion techniques can deliver 
useful information for these more complicated models. 
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APPENDIX: 

The full strong coupling series up to order are 

_ 3a;^ 2 (29 + 20/x) (1443 + 2000/i + 700 

+ 2/i ^ {l + 2fif ~ + (3 + 2/r) ^ {1 + 2 fiY {3 + 2 fif 

4 (51093 + 126936/i + 117230/i2 +47484/i3 + 7056/i^) 

{1+ 2 fif {3 + 2 fif (5 + 2/i) 

^ ^ 2a;2 2 (5 + 2/i) 4 (59 + 68/i + 24/i^ + 4/i^) 

— ujq — 1 -|- 2 -|--—- o — -|-- F - 

l + 2/i (1 + 2/i)^ (1 + 2/i)® (3 + 2/i) 

2 (3313 + 6056/i + 3934/i^ + l 188/i^A 216/i^ +16/i^) 

(l + 2/i)^(3 + 2/i)" 

+2(3578209 + 11359410/i + 14982934/i2 + 10681694/i=^ + 4498120/i^ + 1134808/i^ 

+ 148448/i® - 1888/i^ - 3712/iS - 384/1^)0;^°[(1 + 2/i)®(3 + 2/i)^ (5 + 2/i) (7 + 2/i)]-^ (A2) 

. ^ 6a;2 2 (13 + 2/i) 4 (143 + 112/i + 4/i2 - 4/i3) 

l + 2/i (1 + 2/i)^ (1 + 2/i)^ (3 + 2/i) 

, 2 (-7905 - 11632 /i - 4510 /i^ + 28 /i^ + 104 /i^ - 16 /i^) 

^ (1 + 2/i)^(3 + 2/i)" 

+ (17055678 + 48682964/i + 54776904/i2 + 30230052/i3 + 8173504/i^ + 987088/i® 

+ 157952/i® + 74944/i^ + 14336/i® + 768/i^)a;^°[(l + 2/i)®(3 + 2^if (5 + 2/i) (7 + 2/i)]-^ (A3) 
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FIGURES 


FIG. 1. Finite-size scaling behaviour of the vector energy gap {Mi — 2m)/g as a function of 
I/N"^ for m /(7 = 0, X = 4, where N is the lattice size. Solid boxes mark the finite-lattice data 
points; the dashed line shows a polynomial fit of the form (^). 

FIG. 2. Ground-state energy per site uJo/{xN) as a function of y = Ij^/x for m/g = 0, The 
dashed lines are integrated differential approximants to the series data, while the solid boxes with 
error bars represent finite-lattice estimates. The solid line is a fit to the finite-lattice data in powers 
of y. The exact continuum limit is marked by an open circle. 


FIG. 3. Vector mass gap Mijg as a function of y = \fx for mjg = 0. Dashed lines are series 
approximants, solid boxes are finite-lattice data, and the open circle marks the exact continuum 
limit. The solid line is a linear fit in powers of y to the finite-lattice data. 

FIG. 4. Scalar mass gap M 2 /g as a function of y = for m/g = 0. Notation as Fig. 3. 

Here the solid line represents a second-order polynomial fit to the finite-lattice data over the range 
[0.2, 0.8] in y. 


FIG. 5. The vector energy gap (Mi — 2m)/g as a function of y = l/\/x for m/g = 32. Notation 
as Fig. 4. Here the solid line represents a second-order polynomial fit to the finite-lattice data over 
the range [0-0.5] in y. 


FIG. 6. The scalar energy gap (M 2 — 2m)/g as a function of y = l/\/x for m/g = 32. Notation 
as Fig. 4. Here the solid line represents a polynomial fit to the finite-lattice data over the range 
[0-0.7] in y. 


FIG. 7. Estimates of the bound-state energies Ei/g versus m/g for both the vector and scalar 
states. The open circles are our series estimates, the solid boxes are our finite-lattice estimates, the 
dashed line is the result of Mo and Perry, and the solid lines represent the asymptotic behaviour 
in the limiting case m/g 0, and the non-relativistic limit. 

FIG. 8. Graph of the quantity fo/v/^^ as a function of z = l/u for the vector excited state. 
Notation as in Fig. 4. The solid line represents a polynomial fit in 1/n to the series results. 

FIG. 9. Graph of the quantity fo/u^^^ as a function of z = l/« for the scalar excited state. 
Notation as Fig. 4. The solid line is a linear fit in z = l/n to the series data over the range [0.6-1] 
in z. 
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FIG. 10. Graph of the quantity as a function of z = 1/tt for the vector excited state. 

The solid lines are integrated differential approximants to the series data, and the solid boxes are 
finite-lattice estimates. 
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TABLES 


TABLE I. Series coefficients in x for the series /o, /i, / 2 , and /a for the vector and scalar 


excited states. Nonzero coefficients of x'^ are listed. 


n 

fo 

fi 

A 

fs 



vector excited state 


0 

1.000000000000 

0.000000000000 

0.000000000000 

0.000000000000 

2 

2.000000000000 

-4.000000000000 

8.000000000000 

-1.600000000000x101 

4 

-1.000000000000 xlQi 

5.600000000000x101 

-2.160000000000x102 

7.040000000000x102 

6 

7.866666666667x101 

7,484444444444x102 

4.344296296296x10° 

-1.979753086420x10* 

8 

-7.362222222222x102 

9.942962962963x10° 

-7.742029629630x10'* 

4.547973004115x10° 

10 

7.572929100529 xlO® 

-1.326103683144x10° 

1.296902835617x10° 

-9.337217941993x10° 

12 

-8.273669056689 X10-1 

1.780034473768x10° 

-2.096087983379x101 

1.787437096185x10° 

14 

9.428034196036x105 

-2.405247329390x101 

3.312148486580x10° 

-3.264352340840x10° 

16 

-1.108357531764x101 

3.269962317907x10° 

-5.154702578031x10° 

5.764697147243x101° 

18 

1.334636403738x108 

-4.469545711783x10° 

7.935843765541x101° 

-9.927975671713x1011 

20 

-1.637995781331x10® 

6.137741867815x101° 

-1.211925759935x1012 

1.676878382095x101° 

22 

2.041592824001x101° 

-8.462460202654x1011 

1.839263212769x101° 

-2.788695534036x101* 

24 

-2.577314968112x1011 

1.170804957950x101° 

-2.777427143130x101-1 

4.579160745014x101° 

26 

3.288641969992x1012 

-1.624677659161x101'* 

4.176990920608x101° 

-7.439926786605x101° 

28 

-4.234675197596x101° 

2.260335118290x101° 

-6.260315585321x101° 

1.197974308422x101° 



scalar excited state 


0 

1.000000000000 

0.000000000000 

0.000000000000 

0.000000000000 

2 

6.000000000000 

-1.200000000000x101 

2.400000000000x101 

-4.800000000000x101 

4 

-2.600000000000x101 

1.520000000000x102 

-6.000000000000x102 

1.984000000000x10° 

6 

1.906666666667x102 

-1,884444444444x10° 

1.120829629630x10-1 

-5.195753086420x10* 

8 

-1.756666666667x10° 

2.435066666667x10-1 

-1.932472592593x10° 

1.151850271605x10° 

10 

1.804833650794x10-1 

-3.218263683144x10° 

3.192825551667x10° 

-2.325598645777x101 

12 

-1.979052000756x10° 

4.314452729791x10° 

-5.136573365113x101 

4.420974571815x10° 

14 

2.267367521095x10° 

-5.842838007004x101 

8.115577786905x10° 

-8.058540815871x10° 

16 

-2.681007957489x101 

7.972653477157x10° 

-1.265571251416x101° 

1.424081244509x1011 

18 

3.246557142476x10° 

-1.094270467476x101° 

1.954139217902x1011 

-2.457363530118x1012 

20 

-4.005343740541x10° 

1.508952103990x1011 

-2.994043095570x1012 

4.161067716407x101° 

22 

5.016016326942x101° 

-2.088733926173x1012 

4.558824467106x101° 

-6.938732610121x101* 

24 

-6.359434181405x1011 

2.900477258173x101° 

-6.905937254404x101-* 

1.142452124172x101° 

26 

8.145971243132x1012 

-4.038508440296x101-1 

1.041673344896x101° 

-1.861014699378x1011 


TABLE II. Table of estimates of oo, oi, • • •, 03 for the vector and scalar masses, 
from series data, finite-lattice data, and exact results [P,p^. 

obtained 


ao /2 

ai + 2 

2a2 

403 



Vector state 



Series 

0.56(2) 

1.80(2) 

0.16(4) 

- 0 . 22 ( 6 ) 

Finite lattice 

0.57(1) 

1.78(2) 

0.3(1) 


Exact 

0.5642 

1.781 

0.1908 




Scalar state 



Series 

1.25(15) 

3.2(2) 

-4(2) 


Finite lattice 

1.14(3) 

2.5(10) 

0 ( 2 ) 


Exact 

1.1284 

3.562 

-0.485 
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TABLE III. Estimates of the continuum bound-state energies Eijg, E 2 Ig as functions oimjg. 


The series and finite-lattice estimates obtained in this work are compared with the earlier fi¬ 


nite-lattice estimates of Crewther and Hamer 
and Mo and Perry |T( 


and the light-cone estimates of Eller et al. |15] 


m/g 

series 

this work 

finite-lattice 
this work 

Eg 

ffi 

Efier et al. 

[R 

Perry 




Vector state 



0 

0.56(2) 

0.57(1) 

0.56(1) 



0.125 

0.53(3) 

0.52(2) 

0.54(1) 

0.60 

0.54 

0.25 

0.52(4) 

0.52(2) 

0.52)1) 

0.53 

0.52 

0.5 

0.50(2) 

0.50(2) 

0.50(1) 

0.49 

0.49 

1 

0.46(4) 

0.46(3) 

0.46(1) 

0.44 

0.44 

2 

0.41(4) 

0.41(3) 

0.413(5) 

0.39 

0.39 

4 

0.34(2) 

0.35(2) 

0.358(5) 

0.34 

0.34 

8 

0.30(3) 

0.31(2) 

0.299(5) 

0.28 

0.29 

16 

0.24(3) 

0.25(2) 

0.245(5) 

0.23 

0.24 

32 

0.20(4) 

0.20(2) 

0.197(5) 

0.20 

0.20 




Scalar state 



0 

1.25(15) 

1.14(3) 

1.12(5) 



0.125 

1.35(15) 

1.24(4) 

1.11(5) 

1.41 

1.22 

0.25 

1.30(15) 

1.27(5) 

1.12(5) 

1.31 

1.23 

0.5 

1.25(15) 

1.25(5) 

1.15(5) 

1.23 

1.20 

1 

1.10(15) 

1.14(5) 

1.19(5) 

1.13 

1.12 

2 

1.00(15) 

1.01(5) 

1.10(5) 

0.98 

0.99 

4 

0.90(15) 

0.85(4) 

0.93(5) 

0.84 

0.84 

8 

0.70(15) 

0.73(5) 

0.77(5) 

0.69 

0.70 

16 

0.5(1) 

0.59(5) 

0.62(5) 

0.55 

0.56 

32 

0.4(1) 

0.50(5) 

0.49(5) 

0.46 

0.46 
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TABLE IV. Series coefficients in u for the series /o, /i, / 2 , and /s for the vector and scalar 


excited states. 


n 

/o 

h 

h 

h 



vector excited state 


0 

1.000000000000 

0.000000000000 

0.000000000000 

0.000000000000 

1 

1.000000000000 

-5.000000000000x10-1 

2.500000000000x10-1 

-1.250000000000x10-1 

2 

-5.000000000000x10-1 

-5.000000000000x10-1 

1.125000000000 

-1.250000000000 

3 

2.500000000000x10-1 

5.000000000000x10-1 

6.875000000000x10-1 

-3.437500000000 

4 

-6.250000000000x10-2 

-4.375000000000x10-1 

-6.718750000000x10-1 

-1.375000000000 

5 

-4.687500000000x10-2 

2.500000000000x10-1 

6.484375000000x10-1 

1.496093750000 

6 

7.161458333333x10-2 

-3.906250000000x10-3 

-5.179036458333x10-1 

-1.298828125000 

7 

-3.765190972222 x 10-2 

-1.729600694444x10-1 

2.191840277778x10-1 

9.108886718750x10-1 

8 

-1.051613136574x10-2 

1.898419415509x10-1 

1.635357892072x10-1 

-4.226029007523x10-1 

9 

3.595046055170x10-2 

-6.325276692708x10-2 

-4.165481755763x10-1 

-1.153888466917x10-1 

10 

-2.772084302863x10-2 

-9.454201278373x10-2 

3.571410159515x10-1 

5.813870214140x10-1 

11 

6.427749717514x10-"* 

1.602973539774x10-1 

-1.180466277773X 10-2 

-7.371812834432x10-1 

12 

2.087420768938x10-2 

-9.161565485249x10-2 

-3.598007444620X 10-1 

3.895490578301x10-1 

13 

-2.186375166505x10-2 

-4.596852371627x10-2 

4.496844531498x10-1 

3.420860280769x10-1 

14 

5.564121895836x10-3 

1.366704793607x10-1 

-1.603366455727x10-1 

-9.530272768917x10-1 



scalar excited state 


0 

1.000000000000 

0.000000000000 

0.000000000000 

0.000000000000 

1 

3.000000000000 

-1.500000000000 

7.500000000000x10-1 

-3.750000000000x10-1 

2 

-5.000000000000x10-1 

-2.500000000000 

4.125000000000 

-4.250000000000 

3 

-2.500000000000x10-1 

1.250000000000 

3.562500000000 

-1.187500000000x101 

4 

-6.250000000000x10-2 

8.125000000000x10-1 

-4.046875000000 

-4.187500000000 

5 

4.687500000000x10-2 

1.718750000000x10-1 

-2.632812500000 

1.437109375000x101 

6 

7.161458333333x10-2 

-3.554687500000x10-1 

-9.993489583333x10-2 

8.261718750000 

7 

3.765190972222x10-2 

-4.944661458333x10-1 

2.158637152778 

-2.571695963542 

8 

-1.051613136574x10-2 

-2.319607204861x10-1 

2.582842791522 

-1.177163357205x101 

9 

-3.595046055170x10-2 

1.738077799479x10-1 

8.375302538460x10-1 

-1.176848630552x101 

10 

-2.772084302863x10-2 

3.925574031877x10-1 

-1.691848582201 

-9.678567015095x10-1 

11 

-6.427749717514x10-"* 

2.711228939731x10-1 

-2.874944640017 

1.281942563814x101 

12 

2.087420768938x10-2 

-6.325043120993x10-2 

-1.598536929541 

1.730423290807x101 

13 

2.186375166505x10-2 

-3.223758205015x10-1 

1.178489577243 

6.362423273963 


20 
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